When the Lowest Energy Does Not Induce Native Structures: Parallel Minimization of Multi-Energy Values by Hybridizing Searching Intelligences
نویسندگان
چکیده
BACKGROUND Protein structure prediction (PSP), which is usually modeled as a computational optimization problem, remains one of the biggest challenges in computational biology. PSP encounters two difficult obstacles: the inaccurate energy function problem and the searching problem. Even if the lowest energy has been luckily found by the searching procedure, the correct protein structures are not guaranteed to obtain. RESULTS A general parallel metaheuristic approach is presented to tackle the above two problems. Multi-energy functions are employed to simultaneously guide the parallel searching threads. Searching trajectories are in fact controlled by the parameters of heuristic algorithms. The parallel approach allows the parameters to be perturbed during the searching threads are running in parallel, while each thread is searching the lowest energy value determined by an individual energy function. By hybridizing the intelligences of parallel ant colonies and Monte Carlo Metropolis search, this paper demonstrates an implementation of our parallel approach for PSP. 16 classical instances were tested to show that the parallel approach is competitive for solving PSP problem. CONCLUSIONS This parallel approach combines various sources of both searching intelligences and energy functions, and thus predicts protein conformations with good quality jointly determined by all the parallel searching threads and energy functions. It provides a framework to combine different searching intelligence embedded in heuristic algorithms. It also constructs a container to hybridize different not-so-accurate objective functions which are usually derived from the domain expertise.
منابع مشابه
Protein folding by restrained energy minimization and molecular dynamics.
Native-like folded conformations of bovine pancreatic trypsin inhibitor protein are calculated by searching for conformations with the lowest possible potential energy. Twenty-five random starting structures are subjected to soft-atom restrained energy minimization with respect to both the torsion angles and the atomic Cartesian co-ordinates. The restraints used to limit the search include the ...
متن کاملMulti-Component-Multiphase Flash Calculations for Systems Containing Gas Hydrates by Direct Minimization of Gibbs Free Energy
The Michelsen stability and multiphase flash calculation by direct minimization of Gibbs free energy of the system at constant temperature and pressure, was used for systems containing gas hydrates. The solid hydrate phase was treated as a solid solution. The fugacities of all components of the hydrate phase were calculated as a function of compositions by the rearranged model of van der Wa...
متن کاملCan a physics-based, all-atom potential find a protein's native structure among misfolded structures? I. Large scale AMBER benchmarking
Recent work has shown that physics-based, all-atom energy functions (AMBER, CHARMM, OPLS-AA) and local minimization, when used in scoring, are able to discriminate among native and decoy structures. Yet, there have been only few instances reported of the successful use of physics based potentials in the actual refinement of protein models from a starting conformation to one that ends in structu...
متن کاملContribution of energy values to the analysis of global searching molecular dynamics simulations of transmembrane helical bundles.
Molecular interactions between transmembrane alpha-helices can be explored using global searching molecular dynamics simulations (GSMDS), a method that produces a group of probable low energy structures. We have shown previously that the correct model in various homooligomers is always located at the bottom of one of various possible energy basins. Unfortunately, the correct model is not necess...
متن کاملEnergetics of the native and non-native states of the glycophorin transmembrane helix dimer.
Using an implicit membrane model (IMM1), we examine whether the structure of the transmembrane domain of Glycophorin A (GpA) could be predicted based on energetic considerations alone. The energetics of native GpA shows that van der Waals interactions make the largest contribution to stability. Although specific electrostatic interactions are stabilizing, the overall electrostatic contribution ...
متن کامل